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We consider radially twisted nanotubes in the low-energy approximation where the dynamics 
is governed by a one-dimensional Dirac equation. The mechanical deformation of the nanotubes 
is reflected by the presence of an effective vector potential. We discuss twisted carbon and boron- 
nitride nanotubes where deformations give rise to periodic and non-periodic finite-gap Hamiltonians. 
The intimate relation of these systems with the integrable AKNS hierarchy is exploited in the study 
of their spectral properties as well as in the computation of the (local) density of states. We show 
that a nonlinear hidden supersymmetry generated by local supercharges arises naturally in the 
finite-gap configurations of twisted nanotubes with time-reversal symmetry. The properties of the 
van Hove singularities are encoded in its structure. 



I. INTRODUCTION 

Since their discovery in 1991 [H, carbon nanotubes at- 
tract massive attention of both experimental and theo- 
retical physicists. Their remarkable mechanical and elec- 
tronic properties, extreme mechanical strength 0, elas- 
ticity as well as conductivity [1], make their use in 
the future electronic devices very promising Carbon 
nanotubes are also relevant as a low-dimensional test field 
of fundamental physical theories. For instance, the Klein 
tunneling, a well known phenomenon in relativistic quan- 
tum theory 0] , has not been observed for elementary par- 
ticles so far. However, it was predicted and observed in 
carbon nanotubes where it stays behind the absence of 
back-scattering on impurities [3] . 

The single-wall carbon nanotubes are small cylinders 
with atom-thick shell, that can be created by rolling 
up graphene nanoribbons Despite their structural 
simplicity (as quasi one-dimensional objects), the nan- 
otubes exhibit a remarkable variety of electronic prop- 
erties. They can be either metallic or semiconducting, 
dependently on the orientation of the hexagonal lattice 
in the nanotube. The gap between the valence and con- 
duction band can be regulated b y ei ther external fields @ 
or by mechanical deformations [T^l- This paves the way 
to strain engineering (llj where graphene-based devices 
would be fine-tunned by deformations. 

In the present article, we will consider a class of ex- 
actly solvable models of single-wall carbon nanotubes 
subject to radial twist (axial torsion). The systems will 
be studied in the low-energy regime where the dynam- 
ics of charge carriers is well app roximated by a (1 + 1)- 
dimensional Dirac equation [l2|. In particular, the twist 
configurations under investigation will be described by a 
finite-gap Hamiltonian. 

Besides the analysis of the spectral properties, the 
computation of the local density of states (LDOS) of the 
systems will be addressed. This quantity is measurable 
by means of scanning tunneling microscopy experiments 
[l3| and is very important for the specification of the 
electronic properties of nanostructures. When integrated 
spatially, it provides the density of states (DOS) that re- 



fiects the probability of inserting an electron at given 
energy into the system. 

In general, it is a rather complicated task to analyti- 
cally compute the local density of states. However as it 
was suggested in Ij] and demonstrated later in [1^, the 
computation can be considerably simplified for a broad 
class of Dirac and Bogoliubov de Gennes Hamiltonians 
that are characterized by a finite number of gaps in the 
spectrum. This approach is based on the intrinsic proper- 
ties of the finite-gap systems that arise from the station- 
ary AKNS (Ablowitz-Kaup-Newell-Segur) hierarchies of 
integrable systems (Toj . 

In the next section, we will set up the theoretical frame- 
work in which the twisted nanotubes will be analyzed. 
It will be explained how the Dirac Hamiltonian emerges 
in the description of twisted carbon (and boron-nitride) 
nanotubes. In section Hill the relevant points of the con- 
struction of the AKNS hierarchies are summarized and 
the closed formula for the LDOS is discussed. The sec- 
tion IIVI is devoted to presentation of explicit examples 
of carbon and boron-nitride nanotubes where the LDOS 
and DOS will be computed explicitly. We will show in 
section |V] that the singular points of DOS, so called van 
Hove singularities, are closely related to the nonlinear su- 
persymmetry that arises naturally for the nanotubes with 
finite-gap twists and conserved time-reversal symmetry. 
Last section is left for the discussion and outlook. 



II. THE MODEL 

Carbon atoms have four valence electrons; three of 
them arc tightly binded in the interatomic bonds while 
the fourth one is free and can contribute to the electronic 
properties of the crystal. The properties of the collec- 
tive excitations of these electrons in graphene are well 
described by the tight-binding Hamiltonian (T^]. The 
interactions between the nearest neighbors are assumed 
only, being specified by a constant hopping parameter. 
The spin degree of freedom of electrons can be neglected; 
it is irrelevant in the considered interactions. 

The analysis of the band structure displays the spe- 
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cific feature of graphene: the Fermi surface is formed by 
six discrete points where valence and conduction bands 
meet. They are located at the corners of the hexago- 
nal first Brillouin zone and are called Dirac points. Only 
two of them are inequivalent^ and correspond to different 
electronic states. Let us denote them as K and K' = — K. 
In the close vicinity of the Dirac points, the energy sur- 
face acquires cone-like shape. It suggests that the dis- 
persion relation is linear in this region. Indeed, taking 
k = ±K -I- 5k with (5k ~ and expanding the tight- 
binding Hamiltonian up to the terms linear in (5k, we get 
the stationary equation for the two-dimensional massless 
Dirac particle [12[. The Hamiltonian acquires the same 
form in both valleys of K and K'. In coordinate repre- 
sentation, we have 

h{±K)tp±K = {~i(^2dx + icri9y)V'±K = £'V'±k ■ (1) 

The spinorial degree of freedom in ([l}, the pseudo-spin, 
arises due to the two carbon atoms in the elementary cell; 
the hexagonal lattice can be thought as assembled from 
two triangular lattices. The spin- up or down components 
of the wave functions are non-vanishing only on one of 
the two triangular sublattices. The operator /i(K) acts 
on the spinors (V'ka, V'kb)* while h{K.') on the spinors 
("^K'S, V'k'a)*- Here, the first index denotes the valley, 
the second distinguishes between the sublattices A and 
B and t denotes transposition. 

The formula (H]) was introduced already in 1984 by 
Semenoff and makes the basis for considering the 
condensed matter system as a convenient test field for a 
low-dimensional quantum field theory. Indeed, it makes 
possible to observe phenomena in this condensed matter 
system that are native in high-energy quantum physics, 
see e.g. El, [H. 

A single-wall carbon nanotube is rolled up from a 
straight graphene strip. The actual orientation of the 
hexagonal lattice in the strip is uniquely determined by 
the chiral (circumference) vector Ch which is a linear 
combination of the translation vectors of the lattice [l^l . 
Its length corresponds to the diameter of the nanotube. 
We can fix the coordinates such that y goes in the cir- 
cumference direction. Then the chiral vector gets simple 
form, Ch = (0,C,,). 

The effect of rolling up the strip is reflected by the 
periodic boundary condition imposed on the wave func- 
tions, V'K(a;, y + Gh) = ipK.{x, y). It leads to the quanti- 
zation of the momentum in the circumference direction 
which acquires discrete values ky. In the low-energy ap- 
proximation, only the value of ky that minimizes the en- 
ergy is relevant. The system is then governed by a truly 
one-dimensional Hamiltonian ~i(J2dx +criky. The actual 
value of this fixed ky depends on the character of the 



nanotube. Instead of going into more details that can be 
found for instance in Ref. [l^l, let us notice that ky = 
corresponds to metallic nanotubes as there is no gap in 
the spectrum. When ky ^ 0, there is a small gap in 
the spectrum and the nanotube is semiconducting. For 
purposes of our current analysis, we can suppose that 
the nanotubes are metallic (i.e. the angular momentum 
is vanishing, ky — 0) and are infinitely long. The lat- 
ter approximation is rather reasonable due to the recent 
experiments where ultra long single- wall nanotubes were 
created 
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Up to now, we considered systems where neither exter- 
nal fields nor any strains were present. By deforming the 
crystal mechanically, the interatomic distances in the lat- 
tice are modified. Thus, the hopping parameter ceases to 
be constant and becomes position-dependent. This leads 
to the appearance of gauge fields in the tight-binding 
Hamiltonian. It can be approximated in the low-energy 
limit by Dirac operator with non- vanishing vector poten- 
tial [ia,[2i,[2|. 

We shall consider radial twist of the nanotubes. Let us 
mention that both single- wall and multi-wall twisted car- 
bon nanotubes were prepared experimentally by clamp- 
ing small paddles to the ends of the nanotubes, which 
were then rotated [l^l . 

The radial twist shifts the atoms in the lattice perpen- 
dicularly to the axis, preserving the tubular shape of the 
nanotube. The displacement is reflected by a deforma- 
tion vector which measures the difference between actual 
and equilibrium position of atoms. It can be written in 
our specific case as d = (0, dy[x)Y . We consider the situ- 
ation where the displacement is smooth and small on the 
interatomic distance. Then the interaction does not mix 
the valleys of K and K' and the system can be studied in 
the vicinity of one Dirac point only. Stationary equation 
for low energy Dirac fermions in the K-valley acquires 
the following simple form [ll , [l^ , 



/l(K)(/. = (-2(729^ + S(.t)(Ti)0 = A(/) , 



(2) 



where the vector potential S(a;) refiects the twist. It is re- 
lated with the displacement vector by d = (0, / 'T.{x)dx). 
In this framework, the constant vector potential S(a;) = 
7 > would refiect a linear displacement d = (0,7a;). 

Finally, let us consider the following generalization of 
(PI where a mass term is included. 



h{'K)(l> = {-i(72dx + i;(a;)cri + Mas 



M. (3) 



^ The remaining four Dirac points can be reached by translational 
vectors of the reciprocal lattice, and, hence, do not represent 
distinct electronic states. 



The analogue of ([T|) with the mass term M was pro- 
posed by Semenoff for description of the quasi-particles 
in the boron-nitride crystal in the low-energy approxima- 
tion [T^]. The boron- nitride crystal has the same hexag- 
onal structure as graphene. However, the atoms in the 
elementary cell of the crystal cease to be equivalent. It 
gives rise to the potential term with 0-3 that distinguishes 
between the two triangular sublattices A and B. 

Depcndently on the value of M, we will consider ^ as 
the Hamiltonian of radially twisted, either carbon {M = 
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0) or boron- nitride (M 7^ 0), nanotubes. Boron- nitride 
nanotubes were studied theoretically and observed ex- 
perimentally, see e.g. [1^, [13]. Contrary to the carbon 
nanotubes, they are always semiconducting. 

We shall consider the scenario where the term Y,{x)ai + 
Ma^ in ([3]) belongs to the broad class of the finite gap 
potentials. In the next section, wc will show how the 
peculiar properties of finite-gap systems can be utilized 
for computation of the local density of states. 

III. FINITE-GAP TWISTS AND THE LDOS VIA 
AKNS HIERARCHY 

We review here some properties of the integrable 
ANKS hierarchies associated with the Dirac Hamiltonian 
([3]). To explain more easily the main features, we use the 
following unitarily transformed Hamiltonian, 

h = exp ^(K) cxp (i^) 

- A{xr id, ) ' 

where A = 5](x) -I- iM . This form with diagonal deriva- 
tive term, the Bogoliubov de Gennes type Hamiltonian, 
is used frequently in the analysis of Gross-Neveu and 
Nambu- Jona Lasinio models [l^ , [l^ and will make the 
presentation more coherent with the specialized litera- 
ture di]. 

The vector potential A(a;) = A in Q is called finite- 
gap (or algebro-geometric in the mathematical literature) 
when it solves one of the equations from the stationary 
AKNS hierarchy of the nonlinear differential equations, 
namely AKNSjv- One of the most intriguing properties 
of the Hamiltonian Q with a finite gap potential is man- 
ifested in it s sp ectrum: it consists of a finite number of 
bands [l^ , [23 • The actual number of bands (or gaps) is 
fixed by the AKNSa? equation solved by A. The values of 
band-edge energies of a finite-gap system can be obtained 
in purely algebraic manner, see [l6|. These features are 
intimately related with existence of an integral of motion 
of the Hamiltonian Q . Non-periodic finite-gap systems 
can be obtained as the infinite-period limit of the peri- 
odic ones. In this context, the non-periodic systems are 
known as kink or kink anti-kink models in analogy with 
the soliton solutions in the Korteweg-de Vries (KdV) hi- 
erarchy. 

Another relevant feature of this class of models is that 
they can approximate very well any condensed matter 
systems described, in the low-energy approximation, by 
the Hamiltonian ([2]) with a generic periodic potential. 
The Hamiltonian with a generic periodic potential has 
an infinite number of spectral gaps, the width of which 
decreases rapidly with the increasing absolute value of 
energy. Hence, the spectrum of such system can be fitted 
well by a finite-gap one. 

The stationary AKNS hierarchy of nonlinear differen- 
tial equations can be constructed in terms of a Lax pair 



which consists of the Hamiltonian h and a matrix differ- 
ential operator Sm+i^ defined as 

N+l , f \ 

=^Y.( 'r'^' a ) '^^^^ ' e N . (5) 

The functions fnix) and gn{x) are defined recursively in 
the following manner, 

fn = -^/,:_i+A5„, (6) 

5; = z(A*/„_i-A/,:_0 , (7) 
50 = 1, /-i = 0. (8) 

The functions /„ and g-n depend on A (a;) and its deriva- 
tives and also contains integration constants that appear 
due to the integration of see Ref. [l^ for details. 

The operator ([5]) satisfies the following commutation 
relation for any positive integer N , 

[Sn+iM = ^^^0+') , (9) 

The stationary AKNS hierarchy of nonlinear differential 
equations is then defined in terms of the vanishing com- 
mutator ([HI), 

AKNSat^/at+i^O. (10) 

The Hamiltonian h and the operator S^+i are called 
the Lax pair of the stationary AKNS hierarchy. When a 
function A satisfies the (A^-f- 1)"^ order differential equa- 
tion AKNSjv , all the next equations of the hierarchy with 
greater values than (and with the integration constants 
fixed appropriately^) are immediately solved. 

The operators h and Sn+i satisfy the remarkable al- 
gebraic relation, 

2N+1 

~si+i = n - ' (11) 

n=0 

where A„ are band-edge energies. The operator valued 
function on the right hand side is known as the spectral 
polynomial. The integral Sm+i annihilates all the sin- 
glet eigenstates (/)„ of h corresponding to the band-edge 
energies, [h - A„)0„ = 0, 

^jv+i0« = O, n = 0,l,...,2A+l. (12) 

The local density of states p{x, A) is defined in terms 
of the trace of the Green's function, R{x, A) = G{x, x, A), 

pix,X) ^ lim ImTr^(a;,A), (13) 

71" ImA^o+ 



^ The AKNSjv can be written as a linear combination '^ifl 
where the functions /; are defined like in (|8} but with all the 
integration constant that emerge in Q fixed to zero. 
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where the trace is computed over matrix degrees of free- 
dom. The function R{x, A) is also called diagonal resol- 
vent or Gorkov resolvent. 

The spatial integration of LDOS leads to the formula 
for DOS. In case of periodic quantum systems, the inte- 
gration can be performed over one period L ^, 



DOS = J I p{x,X)dx. 



(14) 



Explicit calculation of the Green's function can be 
quite difficult. Nevertheless, the definition (fT3|) suggests 
that the need of its explicit knowledge might be avoided; 
only the diagonal resolvent is required to find LDOS. This 
fact was utilized in [ij and further developed in [l5| . In- 
deed, an exact form of the diagonal resolvent was found 
for a wide class of Hamiltonians The approach was 
based on the fact that R{x, A) has to satisfy the Dikii- 
Eilenberger equation [soj . 



d ~ 



A 

A*{x) 



-A(x) 
-A 



, R{x; A) (73 



= 0, 
(15) 



where A belongs to the spectrum of h. Additionally, 
i?(x. A) has to satisfy the following requirements. 



R^R', detR 



1 



(16) 



where the latter one fixes the normalization of R. For 
more details on the properties of R and derivation of 
P^ . see e.g. the appendix in [sij . 

Making the following ansatz for the diagonal resolvent 

El, 



N+l 



Rix;X) = ^/3„(A) 



fn-li^) 9n{x) 



(17) 



Making the inverse transformation ([4]) , we can find the 
Lax operator associated with the finite-gap Hamiltonian 
© as 

•Sw+i = cxp (^i^^^ S'tv+i exp (^-i^^^ (19) 

JV+l 

= -2 ^ (ffAT+i-;! + (Talm/w-i + ctiRc/at.;) ■ 

(=0 

The diagonal resolvent for the Hamiltonian h can be ob- 
tained directly from ([3]) , since the trace of an operator is 
invariant with respect to similarity transformations. 



IV. EXACTLY SOLVABLE MODELS OF THE 
TWISTED NANOTUBES 

The periodic systems described by ^ can be classified 
in terms of a quantity which we call average twist. It is 
defined in the following manner. 



max(S) + min(S]) 



(20) 



and corresponds to the value around which the potential 
is oscillating. We will present two and four-gap systems, 
denoted as S(x) — A2 and S(a;) = A4 respectively, where 
the average twist is vanishing. Then we will consider 
two simple cases where it acquires non-zero values. They 
correspond to the one, S(a;) = Ai, and three-gap S(a;) = 
A3 systems. The mass term will be identically zero in all 
these models, M = 0. We will see that the actual value 
of the average twist is in correlation with the qualitative 
spectral properties of these models. 

As the last example, we will consider a non-periodic 
system with a constant mass, M 7^ 0. It will serve for 
illustration of a twisted boron-nitride nanotube. 



and substituting (|17p into ([T5|) . the Dikii-Eilenberger 
equation transforms into the two (mutually conjugated) 
non-linear differentials equations of the form of the 
AKNS hierarchy. The diagonal entries in vanish 
identically due to the recurrence relations ([5]). The re- 
sulting equation can be written as 

N+l N+l 

MX)fn -XJ2 MX)fn-l = , (18) 
n=0 n=0 

which can be solved by fixing properly the constants 
/3„(A), see footnote 2 and Ref. [la]. It can be showed 
that the ansatz ([TT]) fulfills the requirements ((T5)) . 



^ For the non-periodic settings, the spatial integration can be di- 
vergent. 

* The ansatz I I17I I is manifestly hermitian. Additionally, it also 
satisfies the second condition in l|16|l . Indeed, one can check 
directly that (dct/?(x, A))' =0 with the use of I II8I I. 



A. Configurations with zero average twist 

1. Two-gap system 

First, let us consider the system governed by ^ with 
the vector potential 



mk 



I sn mx cn mx 
dn mx 



(21) 



where m is a real parameter and k € (0, 1). This vector 
potential is induced by the deformation specified by the 
following displacement vector, see Fig. [U 



(0, — Indnmx) . 



(22) 



The crystal kink two-gap potential ((2T|) is given in 
terms of doubly periodic Jacobi elliptic functions depend- 
ing on the modular parameter k. It has a real period 
L = 2K{k), where K{k) is the complete elliptic integral 



FIG. 1: Twisting for the 2-gap (upper), 3-gap (middle) and 4- 
gap (lower) configurations of carbon nanotubes. In the system 
without twist, the black line would be straight. 



of the first kind. For the definitions and pro perties of 
the elliptic functions, we recommend Refs. [H, [s^. The 
infinite-period limit (fc — > 1) of (PT|) is called the single 
kink vector potential A2 = tanhx. The potential A2 
vanishes in the limit when the modular parameter goes 
to zero. 

The spectrum of the one-dimensional Dirac Hamilto- 
nian /i(K) has two gaps located symmetrically with re- 
spect to zero. The band-ed ge energ ies are Aq = —A3 = 
— m and Ai — — A2 = — m\/I — fc^. The corresponding 
eigenstates ((/i(K) — A„)(/i„ = 0, n = 0, 1, 2, 3) are 



•0 = — snm.x, 



cn mx \ * 



dn mx 



63 = CT3</'o , (23) 



1 Sn TTIX 

01 = I .^ cnmx, , ) , 02 = o's'/'i • (24) 



dm 



The band-edge energies are non-degenerate, while the en- 
ergies from the interior of the bands are doubly degener- 
ated. 

Using directly the formula ([T7)) for TV = 2, we can find 
the explicit form the diagonal resolvent. Its trace then 
reads 



Tri?2(2:;A) 



A^ + lAl 



^(m2-A2)(A2+TO2(fc2_l)) ' 



(25) 



and the associated density of states acquires the following 
form 

DOS2 = hm Ini = , 

TT ImA^o+ ^/{r^n? - A2)(A2 + m2(fc2 - 1)) 

(26) 

where we have used Eq. and ([14]). Notice that DOS2 
is identically zero when A belongs to the prohibited gaps. 
The function in the argument is purely real for these 
values of A and, thus, the imaginary part is vanishing 
identically, see Fig. [2j 



FIG. 2: Illustration of DOS (ge]) of the Hamiltonian Q with 
T,{x) = A2 and Af = with m = 1.5 and k = 0.7. 



the vector potential 



6mk 



sn mx cn mx dn mx 
1 + fc2 + (5 — 3fc2sn2 m 



(27) 



where 6 = y/l — k'^ + k'^. The crystal kink four-gap po- 
tential ([?7)) is an isospectral deformation of the crystal 
kink potential A4/ = 2mfc2 . Both potentials 



reduce to the single kink A4 = A4' = 2TOtanhma; when 
fc — >■ 1. The associated displacement vector in this case 
takes the following form 



d = (0, - ln(l + fc2 + (5 - 3fc2sn2 mx)) 



(28) 



and is illustrated in Fig. [T] The spectrum of ft.(K) has 
five bands and eight band edge states n = 0,...,7, 
which can be defining with help of and operator V = 



^ -I- A4 as follows ^ 



An 



{h{K) - X„ 



0. 



(29) 



Having in the mind the spectral symmetry A -;->■ —A 
(which is valid for any model ([3]) with M = ), it is 
sufficient to present just the first four eigenstates. They 
arc given by 

^0 ~ m (1 + k'^ — 5 — 3fc^sn^ mx) , Aq = —2m\/S , 

(30) 

-01 = cnxsna;, Ai = —m\/ 2 — fc2 + 2(5 , (31) 
■02 dnxsnx, A2 -r77V2fc2 - 1 + 2(5 , (32) 
-03 = cnxdnx, A3 = —m\/ 2(5 — 1 — fc2 . (33) 



2. Four-gap system 



As the next example, we shall consider the 2K{k)- 
periodic system described by the Hamiltonian 1^ with 



This way to express the eigenfunctions is just the essence of usual 
supersymmetric quantum mechanics applie(i for finite-gap poten- 
tials. To avoid the details here, we refer to 1341. |35||. 
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The remaining band edge states that correspond to pos- 
itive energies can be obtained by multiplying the pre- 
sented ones by 0-3. 

The local density of states can be computed using the 
method described in the preceding section, 



Tri?4(x; A) 



ai 



a2M 



Af 



2A4A'I 



8V(Ao-A2)(A2-A2)(A2-A2)(A2-Ai)' 

(34) 

where the constants are 

ai = 8(A'* - Sm^JA^ -I- WJ), (35) 
a2 = A{6^-5m^d). (36) 

The explicit (analytical) form of the density of states is 
rather cumbersome. In Fig. |3l we present the numerical 
computed DOS of the current four-gap system. 

DOS A4 




FIG. 3: Illustration of DOS of the Hamiltonian Q with 
E(a;) = A4 and M = with m = 1 and k = 0.6. 

Let us note that for the system described by the 
isospectral potential A4' , the resolvent trace 



Tri?4/(x; A) 



m + 772 A|, + rjsAj, 



V(Ao-A2)(A2-A'f)(A2-Ai)(A2-A^) 

(37) 

can be written just in terms of a polynomial in A4' where 
Tji, ri2 and 771 are constant depending of A. 



The constant potential can be regarded as periodic with 
the period being equal to any real number L. We can 
compute the average twist as Ec = 7. 

The spectrum of the system has two bands separated 
by a gap of the width 2I]c. This suggest on the rela- 
tion between the average twist and the central gap in the 
spectrum. 



2. Three-gap system 

Let us test the suggestion in the case of a more compli- 
cated system. Its Hamiltonian ^ has the 2/4r(fc)-periodic 
vector potential 



cn 6 dn 6 



+ fc^ sn 5 sn {x) sn [x -\-b) . 



(39) 



which is called crystal kink-antikink three-gap potential, 
[T5| . The real parameter b G (0,A'(fc)) represents the 
distance between the kink and the antikink. 

The vector potential is induced by the displacement 
d = {Q^F{x)) where F{x) is as follows, 

cn b dn \) 1 

F(x) = Il(ksv?b, am x\k) In (l-fc^ sn^&sn^x) . 

sno 2 

The function n(a;x|(/)) is the incomplete elliptic integral 
of the third kind and am x is the Jacobi amplitude. See 
Fig. [T]for illustration. 

When b ~ K{k), (jSH) is reduced to the two-gap vector 
potential (PT|) . In the infinite period limit, the single 
kink-antikink solution is recovered [3^, limfc_>.i Ai{x) = 
coth6 + tanha; — tanh(.T + b). 

The spectrum of © with ([M)) contains three gaps po- 
sitioned symmetrically with respect to zero. The three 
band-edge states with negative energies are 



(/)o ~ (— sn(x), sn{x + b)f , Aq = 
(/ii = (— cn(x), cn{x + b)Y , Ai = 
4>2 = (-dn(x), dn (a- -!-&))* , A2 



1 

sn6 ' 
dn6 

sn6 ' 
cn6 

sn5 



(40) 



B. Semi-conducting carbon nanotubes via non-zero 
average twist 

1. One-gap system 

The simplest example of a finite-gap system with the 
nonzero average twist is given by the Hamiltonian 
with the constant vector potential S(a;) = 7. The two 
band-edge energies correspond to Ao = — Ai = —7. The 
local density of states can be found in the following form 



Tri?4 = 



A 



(38) 



The positive energy states are obtained as 0n+3 = cr^cjin 
and correspond to the energies An+a = — A„ where n = 
0,1,2. 

The trace of diagonal resolvent can be computed from 
p?)) for = 2 in the following form 



TrR3{x;X) 



v/(A^-A2)(A^-A2)(A^-A2)' 



(41) 



where a = 1 + fc^ -I- 2A^ 



sn-^ b ' 



The actual integration 



of the formula above, needed for analytical form of DOS, 
is rather complicated. We present Fig. 2] of DOS for the 
three-gap case that was obtained numerically. 
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The average twist of associated with the potential ((39 
can be found as 



as 



cnb 
snb' 



(42) 



Checking the corresponding band-edge energies A2 and 
A3 in (|40p . we can see that that the gap between the 
positive and negative energies is exactly of the width 2Sc. 

Comparing the spectra of the presented systems, we 
can see that the non- vanishing average twist (|20|) is pro- 
portional to the magnitude of the central spectral gap in 
the system. In the two and four-gap systems, the average 
twist is vanishing and there is no gap between positive 
and negative energies. These nanotoubes are conducting 
in the sense that infinitesimal excitation is sufficient to 
kick the electrons from valence band to conduction band. 
The systems with the non- vanishing average twist are dif- 
ferent. They have a gap between positive and negative 
energies that are equal to 2I]c and, hence, are semicon- 
ducting. 

DOS A3 





FIG. 4: Illustration of DOS of the Hamiltonian Q 
T,(x) = A3 and M = with k = 0.2 and b = 1.5. 



with 



tpQ = ( tanh I sin 2 ^ j ' ~ "^"^^ ^ 



Ao = -1, (44) 



ipi = ( sech ( sin -.t ) , 



-02 = ( tanh I sin -x j , tan - 



Ai = cos -, 
2' 



A9 = 1. 



(45) 
(46) 



The trace of the diagonal resolvent can be computed in 
the following manner, 



A^ — A cos I — i sin^ | sech^ (sin ^x) 



Vr^A2(A-cos|) 



(47) 



It is worth noticing that in the examples of the car- 
bon nanotubes, the trace of diagonal resolvent (|25p . (|34p . 
(|37l) , and (|4ip could be written in terms of the finite gap 
potential and its derivatives. A similar result is known 
for the Schrodinger systems with Lame potential. The 
trace of the diagonal resolvent corresponded in that case 
to a polynomial of the finite-gap potential [s^l • This is re- 
lated with the fact that the square of Dirac operator with 
the four-gap potential A4' corresponds to an extended 
Schrodinger operator with two-gap Lame potential. 

In current case with the non- vanishing mass term, the 
trace of diagonal resolvent can be written as a function 
of the amplitude of the corresponding complex potential 
A (see (HI)), where |Ap = S(x)2 + Af^. For (gS]), we can 
write 



TtR{x;X) = 



P2n{x,9) 



ViV^~A^n„=i(A-A„) 



Here, 



C. Boron- nitride nanotubes 

As the last example, we shall consider a non-periodic 
system with the nontrivial mass term. We take the po- 
tential term of (jS]) in the following form 



S(x) = A^sin — tanh(sin — a;). 



M 



-NCOS-, (43) 



where N is positive integer and is real parameter. The 
potential is classified as (A^ -I- l)-gap as it solves the cor- 
responding equation of the AKNS. It has A^ -I- 2 singlet 
states in the spectrum. Two of them correspond to the 
energies that form the threshold of the continuum spec- 
trum, the rest is associated with bound states of the sys- 
tem. In the case of A^ = 1, the eigenstates are then given 



/ 9 9 \^ 

P2n{x, 9) = J2 C"(A) I sin' - sech' (sin -x) j (48) 

with Cri(A) being specific constants ^. 



For N = 2 the band edge energies are Aq 



-A4 



-2, Ai 



-A3 



/ 7+cos e 



and A2 = 2 cos ^. The constants in 1481 1 are 



then given in the following form 

co(A) = A"* - 2A^ cos - ^X^icosie) + 7)-l- 



— A 15 cos , 
2 V V2 



361 



ci(A) 



- A^ + A cos ( - j + cos{e) -(- 1 j , C2 = J 
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V. VAN HOVE SINGULARITIES AND THE 
NONLINEAR SUPERSYMMETRY 

The densities of states have a set of singular points 
that are called van Hove singularities. A closer inspec- 
tion of the corresponding formulas ([25]), ([Ml), dSH]), (|4T1) . 
and (|47|) shows that the number as well as the position of 
the van Hove singularities coincide precisely with the sin- 
glet band-edge energies of the finite-gap systems. In this 
section, we will show that this coincidence is reflected by 
a non-linear supersymmetry that underlies the finite-gap 
configurations of the twisted nanotubes. 

Quantum systems in presence of a magnetic field cease 
to be time- reversal invariant. The time-reversal opera- 
tor changes the sign of momentum while it preserves the 
coordinate. It changes the sign of the magnetic field. It 
can be represented by an anti-unitary operator T which 
satisfies T^T = 1, TiT = —i and = —1. The latter 
equality arises due to the half-integer spin of the consid- 
ered particles. 

One can check that the anti-unitary operator (72T (T 
denotes complex conjugation, = 1) does not commute 
with the Hamiltonian ^ due to the symmetry break- 
ing term S(a;) (and M). However, we have to keep in 
mind that these terms arise from the tight-binding model, 
which, despite the deformations of the crystal, is timc- 
revcrsal invariant jl7| . 

The time-reversal symmetry of the system in the low- 
energy regime emerges when dynamics in the both valleys 
corresponding to K and K' is taken into account. The 
total Hamiltonian reads 

« - ( T ,„-"k, ) . («) 

where the energy operators ft.(±K) of the subsystems are 
given as 

h{±K) = -ia^d^ ± S(x)cri ± Mcrg . (50) 

The operator (|49|) acts on the bispinors = 
(V'KA, "^KS, '/'K'B, V'k'a) where we use the notation in- 
troduced in the second section below ([l}. 

The Hamiltonian T-L commutes with the time-reversal 
operator T which is defined in the following manner'', 

[•H,r]=0, r-ai^aaT. (51) 

As the considered system consists of a single fermion, 
the Kramer's theorem applies; all the energy levels of 
(|49)) have to be at least doubly degenerate. In case of pe- 
riodic system, the band structure of ft.(K) is determined 



In I19II . the real spin of electrons in taken into account. There, the 
time reversal operator is defined as (cri ® ai)TcT2, where the last 
Pauli matrix acts on the spin-degree of freedom of the electrons. 
As we do not consider real spin of electrons in our model, we 
have to define T as in ISTI l to keep = -1. 



by 27V -|- 2 non-degenerate band-edge energies A„. In 
the infinite period limit, the operator has iV -I- 2 singlet 
states. As we can see from ([SO]) . the operators /i(±K) 
are unitarily equivalent, ft.(K) = a2h{—'K.)a2. Hence, % 
has the same band structure as /i(K) but the degener- 
acy is doubled as is required by the Kramer's degeneracy 
theorem. 

Degeneracy of energy levels is reflected by a set of 
integrals of motion that are based on the Lax integral 
Sn+i, see In the individual subsystems governed 

by /i(±K), the degeneracy is associated with two diago- 
nal operators, Qo, and Q3, 

Q^^={^T n . y Q3 = r3oQo, (52) 

where t^q = cts (8) 1. The inter- valley (Kramer's) degen- 
eracy is naturally reflected by the operators Qi and Q2, 

Ql=Ti2Qo, S2=r22Qo, (53) 

where Tab = ® (Th, a,b ~ 1,2. All these operators 
commute with the total Hamiltonian, 

[Qa,H]=0. (54) 

By construction, these operators close centrally extended 
Lie algebra so(3), 

[Qo,Qa]=0, [Qa.Qb\='2ieabcQc. a,6,c=l,2,3. 

(55) 

The existence of the operators ([55)1 is a direct conse- 
quence of the time-reversal symmetry of (j49p . Indeed, 
(jSip implies the unitary equivalence of the valley Hamil- 
tonians /i(±K) and enables the construction of antidiag- 
onal operators (|53p . 

The action of the integrals is quite nontrivial and de- 
termined by the properties of the Lax operator Sn+i- It 
can be inferred from ([T^ that all doublet states '^2-deg 
oiH, corresponding to the band-edge energies A„, are all 
annihilated by the integrals of motion Qa, 

Qa'^2-deg=0, = 0,1,2,3. (56) 

Let us denote by the subscript k and k'(=-k) the 
states that are non-vanishing in one valley only, i.e. 
i (1 ± T30) 4'±K = ^±K- We can find mutual eigenstates 
and 5"^, of the Hamiltonian T-l, the valley- index op- 
erator T30, and the integrals Qo and Q3. They satisfy the 
following relations, 

{n - A)*±(K') = 0. (^30 - 1)*| = (r3o + 1)*|. = 

(57) 

and 

Q,*±=±7A*|, Q,*J^, = ±7A*|, i = 0,3. (58) 

The eigenvalues can be determined from the spectral 
polynomial pT|) as 

2N+1 

7A = yP(A)= n (^-^«)'/'- (59) 

n=0 
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Hence, the operators Qo and Q3 act on the basis of ^'±k 
as the multipUcation by -P(A), i.e. as the square-root of 
the operator ^JP{^-C). As mentioned above, the roots of 
the spectral polynomial (|59p coincide with the van Hove 
singularities of the analyzed finite-gap systems. 

The two antidiagonal operators Qi and Q2 switch the 
valley index, 

2l*K(K') = ±7A*K'(K)' Q2*iK = ±«e7A*^K , 

(60) 

where e = ±. 

The behavior of the operators can be conveniently 
expressed in terms of superalgebra with supercharges 
Qi and Q2- Let us fix the valley index operator T30 
as the grading of the superalgebra. We can denote 
•^1(2) = 2i(2) as fermionic and Bn2) = Qo(3) bosonic 
operators. The superalgebra is nonlinear and contains 
another fermionic operators T12 and r22, 

[U,Ba\ = [n,Fa] = 0, {-Fa, J-fe} = 25abP{n), (61) 
[Ba,Fb\^2l5a2e3bcTc2P{n), (62) 
['Ba,T-22] = -2i52a^l, [^0,^12] = 2i52a^l, (63) 
{Fa.n2}^25abBi. (64) 

The spectral polynomial P{T-i) emerges naturally in the 
anticommutator of the supercharges. It underlies non- 
linearity of the superalgebra and manifests the intimate 
relationship of between the algebraic structure and the 
van Hove singularities of the considered models. 

Let us stress that the superalgebra (|5T|) - (|M)) exists for 
any finite-gap configuration of the twisted nanotubes de- 
scribed by T-L as long as the Hamiltonian possesses the 
time-reversal symmetry. 

The choice of the grading operator was not unique. We 
could use either T12 or T22 equally well; both of them ei- 
ther commute or anticommute with the considered oper- 
ators. Notice that T12 corresponds to the unitary compo- 
nent of the time- reversal ([51]) . Choosing any of them as 
the new grading operator, qualitatively the same superal- 
gebra would be obtained. The operators (f52|) -([53 |) would 
be just permuted in the role of bosonic and fermionic 
generators. 

Let us notice that in examples presented in the previ- 
ous section, the single valley Hamiltonians with the vec- 
tor potentials pip and (f27| commute with the operator 
uyJZ where TZ is the parity®, TZxTZ = —x. 

Hence, the corresponding Hamiltonian (|49)) is commut- 
ing with T^^TZ. The latter operator also commutes with 
T30 and r22 whereas it anticommutes with all the op- 
erators Qa, a = 0, ..,3. It means that ts^TZ could be 
regarded as a grading operator of a distinct. A'' = 4 su- 
peralgebra that would be generated by four fermionic op- 
erators ([5^ and ([55)) . The nonlinear superalgebra of Bo- 
goliubov de Gennes Hamiltonians generated by nonlocal 



The three-gap setting with H39| l has the nonlocal integral of 
slightly modified form, see [Tsll . 



supercharges was discussed in the literature, we refer to 
[m for more details, see also [ssj . 

The formulas for LDOS and DOS computed in the 
third section with the use of the formula ([T7| have to 
be multiplied by four to get the correct form for the cor- 
responding twisted nanotubes. Indeed, we have to take 
into account the valley degeneracy that we discussed in 
this section, as well the double degeneracy of energy lev- 
els due to (real) spin— i of the particle that was neglected 
up this moment. 

Finally, let us discuss briefly the settings where an ex- 
ternal magnetic field is present in addition to the twists. 
The magnetic field breaks the time-reversal symmetry. 
When the vector potential A^g (corresponding to the 
magnetic field parallel with the axis of the nanotube) is 
included into the Hamiltonian, we have 

hi±K) = -id^a2 ± Y.{x)ai ± Ma^ + A„,gcri. (65) 

We can see that as long as mass term M is vanishing 
and either magnetic field or twists are switched on (i.e. 
S(a;)Amg = 0), all the energy levels have even degeneracy 
due to the unitary equivalence of /i(K) and /i(K'). The 
situation changes when both Ti{x) and A^g are non-zero. 
In that case, we can still have a finite-gap configuration 
in one valley described by h{K). However, in the second 
valley the finite-gap potential is violated in general by 
the changed sign of A^g. 

Curiously enough, we can get for each valley a dif- 
ferent finite-gap configuration by the fine-tunning of the 
external field. As an example, let us consider the situ- 
ation when the low-energy dynamics in the K-valley is 
described by 

/i(K) = — i9j;(T2 + (coth5 + tanha; — tanh(a; + fe))o'i, (66) 

which is an infinite-period limit of the three-gap system 
([55)) . Let us suppose that the vector potential in ([55)) is 
induced both by radial twist and by external magnetic 
field, where Y,{x) = i (2 coth5 + tanhx — tanh(a; -I- 6)) 
and Amg = \ (tanhx — tanh(x + b)). Then the subsys- 
tem in the K'- valley is described by 

/i(K') = -idxCJ2 + cothfecri, (67) 

which is just the trivial one-gap system. In the current 
setting, the twist associated with Yj{x) is asymptotically 
uniform but gets changed in the localized region where 
the (asymptotically vanishing) external magnetic field is 
nonzero. The spectrum of the corresponding total Hamil- 
tonian % has two singlet discrete energy levels corre- 
sponding to the bound states and two doubly degenerate 
levels ± coth h corresponding to the threshold of the posi- 
tive and negative continuum. The other energy levels are 
four-fold degenerate. Both LDOS and DOS for this sys- 
tem can be computed with use of the formulas described 
in the third section. 
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VI. CONCLUSION AND OUTLOOK 

The one-dimensional Dirac operator with finite-gap 
potential appears in a variety of physically interesting 
models 



141, m 



l9f, I40|. In the present article, we il- 
lustrated how the machinery of the AKNS hierarchy can 
be used in the analysis of the twisted nanotubes in the 
low-energy regime, particularly, for the computation of 
the local density of states. 

We showed that the finite-gap, time-reversal invariant 
configurations possess a hidden nonlinear supersymme- 
try that is associated with the Kramer's degeneracy of 
energy levels. Physics of these systems, namely the pres- 
ence of the two valleys at K and K' and the preserved 
time-reversal symmetry, is responsible of the form of the 
Hamiltonian which consists of two copies of the (uni- 
tarily) equivalent single-valley energy operators. 

The current situation differs from the quantum models 
with bosonized supersymmetry [4l[, where a non-local 



integral of motion was identified as the grading operator. 
Both the Hamiltonian and its integrals of motion 
(|52l) and ((53)) , forming the centrally extended so(3) Lie 
algebra, can be graded by a local operator, e.g. by the 
valley index operator t^q. This framework represents a 
non-trivial example of the hidden supersymmetry in the 
sense that it naturally emerges within the unextcnded, 
physical Hamiltonian (^^1) . 

The operator was interpreted as the effective 
Hamiltonian of the twisted carbon (or boron-nitride) 
nanotubes. Qualitatively the same operator (without the 
mass term). 



h{K) = -i(J2dx + {ky + Ay{x))ai, 



(68) 



describes a massless Dirac particle that moves in a fixed 
direction in a two-dimensional deformed crystal where 
the deformation depends on one coordinate only. Such 
systems with periodic, piece-wise constant vector poten- 
tials induced by either external field or strains were con- 
sidered in ^2 and their transport properties and con- 
ductivity computed. It suggests to analyze the systems 
described by finite-gap potentials in a similar manner. 



The known results [43| on the relation between one- 
dimensional DOS and the phase of the transmission am- 
plitude could be particularly helpful in this context. 

The electron interaction in the nanotubes in presence 
of periodic (both electrostatic and magnetic) fields was 
considered in [4^ . Spectral properties of quasi-particles 
in presence of strain-induce d p eriodic field in planar 
graphene were discussed in |45| . It would be interest- 
ing to consider the finite-gap configurations with electro- 
static potential. In this context, the mapping between 
the systems with magnetic and electrostatic field [i^ 
could provide an interesting way how to extend our re- 
sults. 

Another interesting problem is related to the recent 
observation of the new kind of massless Dirac fermions 
in graphene within periodic fields [47| . One can won- 
der about peculiar characteristics of these states in the 
systems described by finite-gap Hamiltonians. 

Finally, this results can be extended with the use of 
Darboux transformations [1^ which preserve the class 
of finite-gap potentials. Let us mention that the recent 
analysis of the LDOS and Green's functions for generic 
Darboux-transformcd Hamiltonians of the twisted car- 
bon nanotubes psj would be particularly useful in this 
context. 



Acknowledgments 

The authors would like to thank to Gerald Dunne 
and Dmitry V. Kolesnikov for discussions. F.C. was 
supported by the Fondccyt grant 11121651 and by the 
Conicyt grant 79112034 and ACT-91. F.C. wishes to 
thank the kind hospitality of the Nuclear Physics In- 
stitute of the ASCR and the Physics Department of 
the University of Connecticut. The Centro de Estudios 
Cientificos (CECs) is funded by the Chilean Government 
through the Centers of Excellence Base Financing Pro- 
gram of Conicyt. V.J. was supported by GACR Grant 
P203/11/P038 of the Czech Republic. 



[1] S. lijima. Nature 354, 56 (1991). 

[2] Q. Zhao, M. B. Nardelli, and J. Bernholc, Phys. Rev. B 

65, 144105 (2002). 
[3] B. I. Yakobson, C. J. Brabec, and J. Bernholc, Phys. Rev. 

Lett. 76, 2511 (1996). 
[4] S. J. Tans et al, Nature (London) 386, 474 (1997). 
[5] P. G. Collins and P. Avouris, Sci. Amer. 283, 62 (2000); 

P. Avouris, J. Appenzeller, R. Martel, and S. J. Wind, 

Proc. IEEE 91, 1772 (2003); 

M. P. Anantram and F. Leonard, Rep. Prog. Phys. 69, 
507 (2006). 
[6] O. Klein, Z. Phys. 53, 157 (1929). 

[7] T. Ando, T. Nakanishi, and R. Saito, J. Phys. Soc. Jpn. 
67, 2857 (1998); 



P. L. McEuen, M. Bockrath, D. H. Cobden, Y-G. Yoon, 

and S. G. Louie, Phys. Rev. Lett. 83, 5098 (1999); 

T. Ando, H. Suzuura, Physica E 18, 202 (2003); 

V. Jakubsky, L. -M. Nieto and M. S. Plyushchay, Phys. 

Rev. D 83, 047702 (2011). 
[8] O. O. Kit, T. Tallinen, L. Mahadevan, J. Timonen, and 

P. Koskinen, Phys. Rev. B 85, 085428 (2012). 
[9] see the chapter E in [2^ and references therein. 
[10] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 78, 1932 

(1997). 

[11] V. M. Pereira, A. H. Castro Neto, Phys. Rev. Lett. 103, 

046801 (2009). 
[12] G. W. Semenofl, Phys. Rev. Lett. 53, 2449 (1984). 
[13] J. Tersoff, D. R. Hamann, Phys. Rev. B 31, 805 (1985); 



11 



Z. F. Wang, Ruoxi Xiang, Q. W. Shi, Jinlong Yang, Xi- 
aoping Wang, J. G. Hon, and Jie Chen, Phys. Rev. B 74, 
125417 (2006). 

[14] G. Basar, G. V. Dunne, Phys. Rev. Lett. 100, 200404 
(2008); 

G. Basar, G. V. Dunne, Phys. Rev. D 78, 065022 (2008). 
[15] F. Correa, G. V. Dunne and M. S. Plyushchay, Annals 

Phys. 324, 2522 (2009) 
[16] F. Gesztesy, H. Holden, Sohton Equations and Their 

Algebro-Geometric Solutions, Cambridge University 

Press, New York (2003). 
[17] A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 

Novoselov, A. K. Geim, Rev. Mod. Phys. 81, 109 (2009). 
[18] A. Shytov, M. Rudner, Nan Gu, M. Katsnelson, L. Lev- 

itov. Solid State Commun. 149, 1087 (2009). 
[19] V. P. Gusynin, S. G. Sharapov, J. P. Carbotte, Int. J. 

Mod. Phys. B 21, 4611 (2007). 
[20] J. -C. Charlier, X. Blase and S. Roche, Rev. Mod. Phys. 

79, 677 (2007). 
[21] Xueshen Wang et al., Nano. Lett. 9, 3137 (2009). 
[22] D. V. Kolesnikov, V. A. Osipov, Phys. Part. Nucl. 40, 

502 (2009). 

[23] M. A. H. Vozmediano, M. L Katsnelson and F. Guinea, 
Phys. Rept. 496, 109 (2010) and references therein. 

[24] E. Joselevich, ChemPhysChem 7, 1405 (2006). 

[25] V. Jakubsky and M. S. Plyushchay, Phys. Rev. D 85, 
045035 (2012). 

[26] E. J. M. Hamilton et.al.. Science 30, 659 (1993). 

[27] A. Rubio, J. L. Corkill, and M. L. Cohen, Phys. Rev. B 
49, 5081 (1994). 

[28] M. Thies, J. Phys. A 39, 12707 (2006). 

[29] C. De Concini, R. A. Johnson, Ergodic Theory Dynami- 
cal Systems 7, 1 (1987). 

[30] G. Eilenberger, Z. Phys. 214, 195 (1968). 

[31] L Kosztin, S. Kos, M. Stone, and A. J. Legget, Phys. 
Rev. B 58, 9365 (1998). 

[32] M. Abramowitz and L A. Stegun, "Jacobian Elliptic 
Functions and Theta Functions." Ch. 16 in Handbook 
of Mathematical Functions with Formulas, Graphs, and 



Mathematical Tables, New York: Dover, pp. 567-581, 
(1972). 

[33] E. T. Whittaker and G. N. Watson, A course of modern 
analysis (Cambridge Univ. Press, Cambridge, 1980). 

[34] D. J. Fernandez, B. Mielnik, O. Rosas-Ortiz and 
B. F. Samsonov, Phys. Lett. A 294, 168 (2002) 
quant- ph/0302204 . 

[35] F. Correa, V. Jakubsky and M. S. Plyushchay, J. Phys. 
A 41, 485303 (2008) arXiv:0806.1614 [hep-th]]. 

[36] R. Jackiw and C. Rebbi, Phys. Rev. D 13, 3398. (1976) 

[37] E. D. Belokolos et al, Algebro-geometric approach to non- 
linear integrable equations, (Springer, Berlin, 1994). 

[38] F. Correa, V. Jakubsky, L. -M. Nieto and 
M. S. Plyushchay, Ph ys. Rev. Lett. 101, 030403 

(2008) arXiv:0801.1671l [hep-th]]. 

[39] J. Feinberg, Annals Phys. 309, 166 (2004). 
[40] A. Arancibia, J. M. Guilarte and M. S. Plyushchay, 
arXiv: 1210.3666 [math-ph]; 

A. Arancibia and M. S. Plyushchay, Phys. Rev. D 85, 
045018 (2012) 

[41] M. S. Plyushchay, Annals Phys. 245, 339 (1996) 
lhep-th /9601116 ; Int. J. Mod. Phys. A 15, 3679 (2000) 
'hep-th/9903130\ 

[42] L. DellAnna and A. De Martino, Phys. Rev. B 79, 045420 

(2009) ; 

S. Gattenloehner, W. Belzig, M. Titov, Phys. Rev. B 82, 

155417 (2010). 
[43] Y. Avishai, B. Band, Phys. Rev. B 32, 2674 (1985). 
[44] D. S. Novikov, Phys. Rev. B 72, 235428 (2005). 
[45] V. K. Dugaev, M. I. Katsnelson, arXiv: 1206.4526 
[46] Liang Zheng Tan, Ch.-H. Park, and S. G. Louie, Phys. 

Rev. B 81, 195426 (2010). 
[47] Ch.-H. Park, Li Yang, Y-W. Son, M. L. Cohen, and 

S. G. Louie, Phys. Rev. Lett. 101, 126804 (2008); 

L. Brey, H. A. Fertig, Phys. Rev. Lett. 103, 046809 

(2009). 

[48] L.M. Nieto, A. A. Pecheritsin, B.F. Samsonov, Annals 
Phys. 305, 151 (2003). 



